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The paper deals with a Dirichlet spectral problem for a singularly perturbed second 
order elliptic operator with rapidly oscillating locally periodic coefficients. We study the 
limit behaviour of the first eigenpair (ground state) of this problem. The main tool in 
deriving the limit (effective) problem is the viscosity solutions technique for Hamilton- 
Jacobi equations. The effective problem need not have a unique solution. We study the 
non-uniqueness issue in a particular case of zero potential and construct the higher order 
term of the ground state asymptotics. 

1 Introduction 

Given a singularly perturbed elliptic operator of the form 

C £ u = e 2 a ll (x,x/e a ) - n + eb> (x, x/e a ) — + c{x,x/e a )u (1.1) 

OXiOXj OXj 

with a small parameter e > 0, we consider a Dirichlet spectral problem 

C £ u = Xu, u = on dfl 
°Mathematics Subject Classification 2000: 
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stated in a smooth bounded domain Q C M. N . We assume that the coefficients a^(x, y), W{x, y) 
and c(x,y) are sufficiently regular functions periodic in y variable, and that a^(x,y) satisfy 
the uniform ellipticity condition. Finally, a > is a fixed positive parameter. Let us remark 
that in the underlying convection-diffusion model e represents characteristic ratio between the 
diffusion and convection coefficients, while e a refers to the microstructure period. 

As well known, the operator C £ has a discrete spectrum, and the first eigenvalue A e (the 
eigenvalue with the maximal real part) is real and simple; the corresponding eigenfunction u £ 
can be chosen to satisfy u £ > in Q. The goal of this work is to study the asymptotic behavior 
of X £ and u £ as e — > 0. 

The first eigenpair (ground state) of ( 11. lft plays a crucial role when studying the large 
time behavior of solutions to the corresponding parabolic initial boundary problem. The first 
eigenvalue characterizes an exponential growth or decay of a typical solution, as t — > oo, while 
the corresponding eigenfunction describes the limit profile of a normalized solution. 

Also, since in a typical case the first eigenfunction shows a singular behavior, as e — > 0, in 
many applications it is important to know the set of concentration points of u £ , the so-called 
hot spots. This concentration set might consist of one point, or finite number of points, or a 
surface of positive codimension, or it might have more complicated structure. An interesting 
discussion on hot spots can be found in |38| . 

Boundary value problems for singularly perturbed elliptic operators have been widely stud- 
ied in the existing literature. An important contribution to this topic has been done in the 
classical work [42J that deals with singular perturbed operators with smooth non-oscillating 
coefficients under the assumption that for e = the problem remains (in a certain sense) 
well-posed. 

The Dirichlet problem for a convection-diffusion operator with a small diffusion and with a 
convection directed outward at the domain boundary was studied for the first time in |19j . The 
approach developed in that work relies on large deviation results for trajectories of a diffusion 
process being a solution of the corresponding stochastic differential equation. 

The probabilistic interpretation of solutions and the aforementioned large deviation principle 
have also been used in j23], [21], [15], where the first eigenvalue is studied for a second order 
elliptic operator being a singular perturbation of a first order operator. 

There are two natural approaches that can be used for studying the logarithmic asymptotics 
of the principal eigenfunction of a second order singularly perturbed operator. One of them 
relies on the above mentioned large deviation results for diffusion processes with a small diffusion 
coefficients. This method was used in [35J for studying operators with smooth coefficients on a 
compact Riemannian manifolds. 

We follow yet another (deterministic) approach based on the viscosity solution techniques 
for nonlinear PDEs. In the context of linear singularly perturbed equations, these techniques 
were originally developed in [18] and followed by [5], |21| . |34j . [9], [10] and other works (see 
also a review in [4]). Since u £ > in f2, we can represent u £ as u £ (x) = e~ Ws ^' e to find that 
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W £ satisfies 



sa l3 (x, x/e a ) 



d 2 W £ 
dxidxj 



+ H(VW £ ,x,x/e a ) = \ £ 



(1.2) 



with H(p,x,y) = a l i(x,y)piPj — ty(x,y)pj + c(x,y), and the Dirichlet boundary condition for 
u £ yields W £ = +00 on dQ. Using perturbed test functions we pass to the limit in (11.21) and 
get the limit Hamilton- Jacobi equation of the form 



with an effective Hamiltonian H(p,x) whose definition depends on whether a > 1, a = 1 or 
< a < 1. We show that in the limit e — > the boundary condition W £ = +00 on dQ yields 



The latter condition is known [40J, [12J as the state constraint boundary condition. Both 
equation (jl.3p and boundary condition f ll.4p are understood in viscosity sense. 

Equations of type (jl.2p have been extensively studied in the existing literature. One can 
find a short review of state of the art in |29j, [25J and in more recent works [8], pQ, see also 
references therein. 

Earlier, singularly perturbed KPP-type reaction-diffusion equations were studied in |32j 
where, in particular, equations with rapidly oscillating coefficients were considered. It was 
shown that the classical Huygens principle might fail to work in this case. 

In the present work, deriving the effective problem f ll.3p - fll.4p relies on the idea of per- 
turbed test functions originally proposed in [16J. We strongly believe that with the help of 
the techniques developed recently in |29J, [26], [30], [1] this result can be extended to a more 
general almost periodic setting as well as random stationary ergodic setting. In other words, 
the periodicity assumption can be replaced with the assumption that the coefficients in ( 11.11) 
are almost periodic or random statistically homogeneous and ergodic with respect to the fast 
variable, at least in the case a = 1. The case a^l looks more difficult and might require some 
extra assumptions. We refer to [JT], [39J , [7J, [29J, [26J, [8], [T] for (far not complete list of) 
various results on almost periodic and random homogenization of nonlinear PDEs. However, 
the essential novelty of this work comes in the (logically) second part of the paper devoted to 
the improved ground state asymptotics and resolving the non-uniqueness issue for ( ll.3p - fll.4p . 
The generalization of this part to non-periodic settings is an open problem. 

Problem fll.3p - fll.4p is known as ergodic or additive eigenvalue problem. Its solvability was 
first proved in [3 lj in periodic setting, more recent results are contained, e.g., in [22j as well 
as in [H], [2] where stationary ergodic Hamiltonians were considered. There exists the unique 
additive eigenvalue A of ( ll.3p - (ll.4p while the eigenfunction W need not be unique even up 
to an additive constant. This non-uniqueness issue is intimately related to the structure of 
the so-called Aubry set of effective Hamiltonian which play the role of a hidden boundary for 
(ll.3p -f jl~4l . Loosely speaking the non- uniqueness in f ll.3p -( !T~4"P appears when the Aubry set is 



H(VW(x),x) = X in ft. 



(1.3) 



H(VW(x),x) > A 



on dQ. 



(1.4) 
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disconnected. By contrast, for every e > the eigenfunction u £ is unique up to a normalization, 
and it is natural to find true limit of W £ = —elogu e among solutions of f ll.3p -( !T74|) . This 
challenging problem is addressed in a particular case of (11. ip with c(x, y) = 0, a > 1. Following 
[36J we introduce the effective drift (convection) and assume that it has a finite number of 
hyperbolic fixed points in Q, and that the Aubry set of the effective Hamiltonian coincides with 
this finite collection of points. It follows from our results that in this case X £ tends to zero as 
e — > 0. We show that X £ /e has a finite limit that can be determined in terms of eigenvalues 
of Ornstein-Uhlenbeck operators in M. N obtained via local analysis of (II .ip at the scale y/e in 
the vicinity of aforementioned fixed points. This, in turn, enables fine selection of the additive 
eigenfunction corresponding to lim e _j.o We- 

2 Main results 

We begin with standing hypotheses which are assumed to hold throughout this paper. We 
assume that Q is connected and has C 2 boundary dfl; the coefficients a u (x, y), V(x, y), c(x, y) e 
C 1 (il x WL N ) are ^-periodic in y functions, where Y = (0,1)^. The matrix {a l ^) i ,- = jjy is 
uniformly positive definite: 

^'(x,2/)OG>HC| 2 >0 VC^O, (2.1) 

and, without loss of generality, we can assume the symmetry = a^ 1 . 

The first eigenfunction u £ of the operator (II. ip can be normalized to satisfy 

1 = maxw e {u £ > in fl), (2-2) 

then its scaled logarithmic transformation 

W £ := -e\ogu £ 

is a nonnegative function vanishing at the points of maxima of u £ . 
The asymptotic behavior of A £ and W £ is described in 

Theorem 1. The eigenvalues X £ converge as e — >• to the limit X, which is the unique real 
number for which problem ( ti.gj) . fll.^| ) has a continuous viscosity solution. The functions W £ 
converge (up to extracting a subsequence) to a limit W uniformly on compacts in Q, and every 
limit function W is a viscosity solution of ( ti.gj) . ( TTjU . 

The effective Hamiltonian H(p,x) in U.3\) is given by the following formulas, depending on 
the parameter a. 

(i) If a > 1 then 

H(p, x) = J H{p, x, y)${y) dy (2.3) 
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where 

H(p,x,y) = a t3 (x,y)p i p j - b j (x,y)p j + c(x,y), 

and$(y) is the unique Y— periodic solution of the equation 9 ®q (a' J (x, y)®) = normalized by 
f Y 0(y)dy = l. 

(ii) If a = 1 then H(p,x) is the first eigenvalue (eigenvalue with the maximal real part) of the 
problem 

d 2 $ d-d — 

a tJ 0, y) dy . dy . + ( x > V) ~ 2qV ( x ' yw) 'Qy~. +H ^ ) ' x > y ^ = H(j) > ^ 2 4 \ 

d{y) is Y -periodic. 

According to the Krein-Rutman theorem H(p,x) is real. 

(Hi) I/O < a < 1 then H(p,x) is the unique number such that the problem 

H(p + V$(y),x,y) = H(p,x) (2.5) 

has a Y— periodic viscosity solution "&(y); here p € M N and x € Q are parameters. 

We note that the effective Hamiltonian H(p,x) is continuous on Mr x Q, convex in p 
and coercive, moreover H(p,x) > mi\p\ 2 — C, mi > 0. The viscosity solutions theory for such 
Hamiltonians is well established. Following |22| (see also [33| ) we present various representation 
formulas for the solutions of problem fll.3p ~ fll.4p . 

Let us rewrite problem fll.3p ~ fll.4J) in the form 

H(VW(x),x)<\ in ft (2.6) 

H(yW(x),x)>\ infi, (2.7) 

i.e. (12. 6p requires that W is a viscosity subsolution in Q while (12. 7p means that W is a viscosity 
supersolution in Q. Then the number Xjj (additive eigenvalue) for which (jl.3p ~ fll.4p has a 
solution is given by 

XjY = inf{A; (125]) has a solution W E C(H)}. (2.8) 

It can also be expressed in terms of action minimization, 

'■t 



1 f- 
Xjj = — lim - inf / L(f), rf) dr, 

i^oo t J 



where the infimum is taken over absolutely continuous curves 77 : [0,t] — > Q, and L(v,x) is the 
Legendre transform of H(p,x), 

L(v, x) = max{f • p — H(p, x)}. 
Let us define now the distance function 
%_ A _(x, y) = sup{W(x) - W(y); W e C(H) is a solution of AM} for A = (2.9) 
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It is known (see, e.g., [22]) that djj_ x _(x,x) = 0, djj_ x _(x,y) is Lipschitz continuous, 
d-jj_ x _(x,y) < djj_ x _(x,z) + djj_ x _(z,y). Besides, for every y G ft the function djj_ Xo (x,y) is 
a solution of ( 12. 6 j) for A = Ajj and H(SJ x d-jj_ x _(x, y), x) > A# in Q \ {y}. The number A^ is 
such that the Aubry set Af7_ A _, 

*<4#_ A _ = {?/ G H; djj_ x _{x,y) is a solution of (l2.7jl for A = A^}, (2-10) 

is nonempty. Note also that the distance function d-jj_ x _(x, y) admits the representation 

dH_ x _(x,y) = inf{^ (1(71,7]) + A^dr, 77(0) = y, V (t) = x,t > o}, (2.11) 
and the Aubry set can be characterized by 

y e A-jj_ x sup inf { f (1(7), 7,) + X-jj) dr, 77(0) = 7)(t) =y,t>s} = 0. (2.12) 

The infimum in (12.111) and (I2.12p is taken over absolutely continuous curves rj : [0,t] — > Q. 

According to the definition of djj_ x _(x,y) every solution W of ( I1.3p -( 1T7^|) satisfies W(x) — 
W(y) < djj_ x _(x, y); this inequality holds, in particular, for all x, y G Ajj_ x _. Conversely, given 
a function g(x) on Ajj_ x _ which satisfies the compatibility condition g(x) — g(y) < d-jj_ x _(x, y) 
Wx, y G Ajj_ x _ then 

W(x) = min{ci 77 _ ATF (x, y) + g(y); y G Ajj_ x _} (2.13) 



is the unique solution of Ql.3p - fll.4p for A = Xj[ such that W(x) = g(x) on Ajj_ x _. In Ap- 
pendix [A] we show the following simple uniqueness criterion for problem ( ll.3p -( TT~4l) : a so- 
lution W (for A = A^) is unique up to an additive constant if and only if Sjj_ x _(x, y) = 
Vx, y G Ajj_ x _, where S-g_ x _(x, y) denotes the symmetrized distance, Sjj_ x _(x, y) = 
dH_ x _(x, y) + d-jj_ x _(y, x). 

The interesting issue of non-uniqueness in the limit (homogenized) problem can be resolved 
by studying next terms in the asymptotic expansion of A e . This question is rather complicated, 
and we mainly focus in this work on a particular case when c(x,y) = and a = 1, so that 
operator (jl.ip takes the form 

C £ u = e 2 aV(x,x/e)—-^- + eV(x,x/e)-^. (2.14) 

OXiOXj OXj 

Moreover we assume that A^- = and that the corresponding Aubry set Ajj has a special 
structure. The analogous result for a > 1 is established in Section [HI 

For a > 1 the effective Hamiltonian H(p, x) is strictly convex in p, i.e. (^ d p 2 g p H(p, %)j 

is positive definite for all p G ~R N and x G f2, see [11], or |13| for a = 1, while for a > 1 the 
Hamiltonian H(p, x) is a quadratic function in p. Note also that if c(x, y) = then H(0, x) = 0. 
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Therefore, the Lagrangian L(v,x) is strictly convex and L(v,x) = max{p • v — H(p,x)} > 
—H(0,x) = 0, in the case c(x,y) = 0, a = 1 we are interested in. Thus we have 



L(v,x) > 0, and L(v,x) = •<=>■ Vj = (0, x). 



On the other hand direct calculations show that 



OH 
o — 1 



— (0,x) = b J (x):= / V(x,y)9*(x,y)dy, (2.15) 
dpj J Y 

the functions V (x) being components of the so-called effective drift b(x) defined by the right 
hand side of (12.15)) via the ^-periodic solution 9* of 

_g_( o « (XlB)ff .)_^(j( (l , B)ff .) = (2.16) 

normalized by j Y 6*dy = 1. (Note that 9* > and it is a C 2 function.) Thus the Lagrangian 
L(v, x) can be represented in the form 

L(v, x) = k 2_,( v j + b J (x)) 2 + L(v, x), where < L(v, x) < k 2^{ v j + (x)) 2 , < k < k. 

This implies, in view of (I2.12p . that the Aubry set A-jj of the Hamiltonian H coincides with 
that of the Hamiltonian ^p 2 — b (x)pj whose corresponding Lagrangian is | XX u j + b\x)) 2 . 
In particular, the additive eigenvalue is zero if and only if there is an orbit r\ : R — >■ Q, 
f] = —b(rj). We moreover assume that 

AhjLQ and Ajj C Q, _ 
Ajj is a finite set of hyperbolic fixed points £ of the ODE x = —b(x). 

Under this assumption we are able to study the leading (of order e) term of the asymptotic 
expansion of A e . This in turn allows us to establish a sufficient condition for selecting the unique 
limit of functions W £ among solutions of the homogenized problem ( II. 3p . ( II. 4p . 

Theorem 2. Let a = 1 and c(x,y) = 0. Then, under condition A2.17\) we have 

\ £ = ea + o{e) ) where a — max{<7(£); £ 6 Affh (2-18) 

and cr(£) is the sum of negative real parts of the eigenvalues of the matrix 

corresponding to the linearized effective drift at £ (since every fixed point £ is assumed to be 
hyperbolic —B(£) has no eigenvalues with zero real part). Moreover, if the maximum in $2.18\) 
is attained at exactly one £ = £ then 
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(i) the scaled logarithmic transformations W £ = —elogu £ of eigenfunctions u £ (normalized by 
\2.2\l ) converge uniformly on compacts in Q to W(x) = d-ff(x, I;), i.e. W is the maximal viscosity 
solution ofH{VW{x),x) = in Q, H(VW{x),x) > on dQ, such that W@ = 0; 

(ii) w £ (£ + y/ez) — > u(z) in C(K) and weakly in H X {K) for every compact K , and the limit u 
is the unique positive eigenfunction of the Ornstain- Uhlenbeck operator, 

Q ij ^~ + ZiB ij jr- = ou tn R N , (2.19) 

OZiOZj OZj 

normalized by u(0) = 1 and satisfying the following condition, u(z)e^ u * sZ ^ 2 ~ u ^ z ^ 2 is bounded 
on ~R N for some /i > and every v > 0. The coefficients in Ii2.19\) are given by B % i = B l i(l;), 
= \ qp Q p . (0; 0/ n s and II U denote spectral projectors on the invariant subspaces of the 
matrix B corresponding to the eigenvalues with positive and negative real parts (stable and 
unstable subspaces of the system Zi = —B^Zj). 

Remark 3. Condition Ii2.17\) is satisfied, in particular, when the vector field b(x,y) is a C 1 - 
small perturbation of a gradient field VP(x) with C 2 potential P(x) having the following prop- 
erties: the set {x G Q; VP(x) = 0} is formed by a finite collection of points in Q and the 
Hessian matrix (^ 9 ^q x . P{x)j a t every such a point is nonsingular (see Appendix W\) . 

Condition Ii2.17\) is satisfied if and only if the vector field b possesses the following properties: 

• b has a finite number of fixed points in Q, say . . . £ n . All of them are hyperbolic, and none 
of them is situated on dQ. 

• Vy G Q, either sup{t < : x y (t) ^. Q} > —oo, or lim x y (t) = £ J for some j G {1, . . . ,n}, 

t— >— oo 

where x y is a solution of the ODE x y = —b(x y ), x y (0) = y. 

• there is no any closed path ^ 1 ,^ 2 , . . . , £ J * = £ jl with k > 2 such that for any two consecutive 
points £ Js and £ Js+1 there is a solution of the equation x = —b(x) with lim x(t) = £ Js and 

t — >— oo 

lim x(t) = Note that £ J1 might coincide with £ J2 . 

t— 5- + 0O 

Remark 4. It is not hard to show that under condition Ii2.17\) we have S-ff(£, £') > V£, £' G 
Ajj, £ 7^ This means that problem U.3\) . does have many solutions unless Ajj is a 

single point. 

Note that condition (I2.17P of Theorem [2] assumes, in particular, that all Co>(and a)-limit 
points of the ODE x = —b(x) are fixed points. Another important case, when the ODE 
x = —b(x) has limit cycles in Q (which is also the case of general position), will be considered 
in a separate paper. 



3 Singularly perturbed operators in the periodic setting 

Consider the spectral problem for singularly perturbed elliptic operators of the form 

C<f er) u = e 2 a ij (x)^-^- + eb>(x)p- + c{x)u, (3.1) 

OXiOXj OXj 
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with F-periodic coefficients a l \ & , c G C 1 (M Ar ), u also being ^-periodic. We assume the 
uniform ellipticity condition a^(x)QCj — m IC| 2 > V£ G Mr \ {0} and the symmetry a u = a Jl . 
Similarly to the case of the Dirichlet boundary condition, the first eigenvalue fi £ of £^ per ^ 
(eigenvalue with the maximal real part) is a real and simple eigenvalue, the corresponding 
eigenfunction u £ can be chosen to satisfy < u £ (x) < maxti e = 1. The asymptotic behavior of 
fi £ as e — » and u £ was studied in [35] using a combination of large deviation and variational 
techniques. We recover hereafter the results of [35] by means of vanishing viscosity approach 
and establish as a bi-product some bounds for derivatives of functions W £ (x) = — elogu £ (x) 
that are essential in the proof of Theorem [TJ 

First we derive the a priori bounds for the eigenvalues. 

Lemma 5. For every e > the eigenvalue jj, £ of Le satisfies the inequalities 

minc(x) < fi £ < maxc(x). (3.2) 
Proof. Let x' be a maximum point of u £ , we have 

V« E (x') = 0, ^ a ^{x')—^-{x') < 0, 

OXiOXj 

therefore c(x')u £ (x') > fi E u £ (x'), i.e. [i £ < maxc(x). Similarly, if x" is a minimum point of u £ 
then [i £ u £ (x") > c(x")u £ (x") and therefore fi £ > minc(x). □ 

Since u £ = e- w ^( x >/ £ we have 

d 2 Wr 8W F dW £ dW £ 

OXiOXj OXi OXj OXj 

The bounds for the first and second derivatives of W £ (x) are obtained in the following 

Lemma 6. There is a constant C , independent of e, such that 

max |W £ | < C, max \d 2 W £ /dx i dx j \ < C/e. (3.4) 

Proof. The proof of the first bound in (13. 4p is borrowed from [18j. Let D\(x) := \VW £ (x)\ 2 
and L*2(x) := Yl \d 2 W £ (x)/dxidxj\ 2 . From (13. 3p in conjunction with (13. 2p we get mDi < 
C(eD^ 2 + D\^ 2 + 1), this in turn implies that 

D x < C(sDl /2 + 1). (3.5) 

Assume that D\ attains its maximum at a point x', then we have VDi(x') = and 
a ij (x')-fjr(x') < or 

^ V)?V)=0 (3.6) 



dxidxk dx 

and 



: y a u ^ d 2 w £ ^ £ y aij d*w £ dw £ &tx , 

' dxidxk dxjdxk dxSxSx^ dxu 

k 



9 



In order to bound the right hand side of (13. 7p we take derivatives of (13. 3p . this yields 

_ ca H 93 Ws - c dalj 92 We _ 2a ij Q2We dWe + ti ^ — + ^L^l _ dc_ 
dxidxjdxk dxk dxidxj dxidxk dxj dxidx^ dx k dx{ dxk 

Then we multiply (13.81) by dW £ /dxk, sum up the equations in k and insert the result into (13.71) 
to obtain 

fw £ d 2 \ 

' dxjdxk dxjdxk 



smD 2 {x') < < c(eD\' 2 {x')Dl'\x') + D 1 (x') + D\'\x')). 



Next we use (13. 5p to get that /^(^O < C/e, and exploiting once more (I3.5P we obtain the first 
bound in ([S3]). 

To show the second bound in (13. 4p we use the following interpolation inequality 

||Vn|||oo < C(\\a lj d 2 u/dx i dxj\\ L °° + \\u\\l°°)\\u\\l°°, (3.9) 

which holds for every F-periodic u with a constant C independent of u. The proof of this 
inequality follows the lines of one in the Appendix of |6] (here it is important that the coefficients 
a y are Lipschitz continuous). We apply (13.91) to (13.81) to obtain 

W&We/dxtdXHtf^ < l^Wc/dx^XjW^ + 1) V/, k, (3.10) 

here we have also used the first bound in (13.41) . From (I3.10p one easily derives the second bound 
in (JS3D- □ 

It follows from Lemma |5] that fi £ — > fi, up to extracting a subsequence. Due to Lemma [6] 
the family of functions W e (x) is equicontinuous, moreover min W e (x) = therefore passing to a 
further subsequence (if necessary) we have W £ (x) — > W(x) uniformly. The standard arguments 
show that the pair \x and W satisfies the equation 

dW dW dW 

a "^wr h,(x) ^ +c(x)= " (311) 

in the viscosity sense. The number /x for which (I3.1ip has a periodic viscosity solution is unique 
(see [3T].|17|). therefore the entire sequence /i £ converges to \i as e — > 0. 



4 A priori bounds 

In this section we show that the eigenvalues X e of ( II. ip are uniformly bounded and the functions 
W e (given by (T5])) uniformly converge on compacts in Q as e — > 0, up to extracting a subsequence. 
We also prove an auxiliary result on the behavior of the minimum points of W £ — <fi (where 
is an arbitrary C 2 function) which is important in the subsequent analysis. 

Because of the Dirichlet boundary condition on the boundary dQ and fast oscillations of 
the coefficients the arguments here are more involved than those in the periodic case. 
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Lemma 7. There is a constant A independent of e and such that 



- A < \ e < supc(x,y). (4.1) 

Proof. The proof of the upper bound follows by the maximum principle as in Lemma [5j 

To derive a lower bound for A e we construct a function v £ and choose a number A > such 
that v £ = on <9ft, and 

C £ v £ -Xv £ >0 in ft (4.2) 

for every A < —A, < e < 1. There is a function W G C 2 (ft) satisfying the following conditions, 
W > in ft and W = on <9ft, |VW| > 1 in a neighborhood of <9ft. Set u e (z) := e KW( - x ^ £ - 1, 
where k is a positive parameter to be chosen later. We assume that —A < min c(x,y) so that 
A < min c(x,y). Then we have 

C e v e - \v £ > (ttik 2 - k(M 1 + eM 2 ) + (c(x, x/e a ) - X))e KWix)/£ - (c(x, x/e a ) - A) > in ft' 

when k > «i := (Mi + M 2 )/m. Here ft = {x G ft; |VW| > 1}, Mi = max\b i {x,y)^(x)\, 
M 2 = max|a u (x, y) ®^ x , (x) | . On the other hand 5 := inf{W(z); x G ft \ ft'} > therefore 
e nW(x)/e > 2 in ft \ ft', when k > k 2 := (log 2) / 5. Assuming additionally that min c(x, y) — A > 
2«(Mi + M 2 ), we have 

C £ v £ -Xv £ > {-K{M 1 +sM 2 ) + {c{x,x/e a )-\))exp{K-f/s)-{c{x,x/e a )-\) > in ft \ ft'. 

Thus setting k := max{/«i, ^2} and A := 2k(Mi + M2) — mmc(x,y) we get f !4 .2D . 

Now note that A £ is also the first eigenvalue (the eigenvalue with the maximal real part) of 
the adjoint operator C* £ u = e 2 Q ® 2 dx {a^u) —e-£r(b l u) + cu, and the corresponding eigenfunction 
u* can be chosen positive in ft. Therefore, if A e < —A then (C £ v £ — X £ v £ )u* > in ft that 
contradicts the Fredholm theorem. □ 

The following two results show that, up to extracting a subsequence, functions W £ converge 
uniformly on compacts in ft. For brevity introduce the notation 

d(x) = dist(x, 9ft). 

Lemma 8. For every k > there is a constant C K , independent of e, such that 

\W £ (x) - W £ (z)\ < C K {\x - z\ + e) when x,z G ft, mm{d(x) , d(z)} > ke. (4.3) 

Proof. Let Xo G ft, d(xo) > ke. Changing the variables x 1— > Xq + ey in operator (II. ip we find 
that v £ = u £ (xq + ey) satisfies the equation 

d 2 v ■ dv 

a7(y) Qy.Qy. + V ^'0y L . + ( C£ ^ ' X ^ VS = When l ^ l < K - 
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Note that a l J , W £ and c £ are uniformly bounded and a l J( l (i > m\(\ 2 > \/( ^ 0. Therefore 
by Harnack's inequality (see, e.g., [27]) we have v £ (y) < C' K v E (y') when \y\, \y'\ < re/2, where 
C' K > is independent of Xq and e. Thus 

| W e (x) - W e (z)\ < (\ogC' K )e, if \x-z\< ne/2 and min{d(x), d(z)} > ke. (4.4) 

Since dfl is C 2 -smooth and Q is connected, every two points x, z G Q such that d(x) > re£, 
> re£ can be connected by a chain of segments i = l,n £ with — Xj| < ree/2, 

dist(xj,9fi) > res and n £ < C"(\x — 2|/e + 1). (We assume that e is sufficiently small.) Then 
iterating (IP)) we obtain P~3)) . □ 

Lemma 9. Let <f)(x) G C 2 (f2) £/ien every maximum point x e of the function u £ (x)e^ x ^ e satisfies 
d(x £ ) > re£ some re > independent of e. 

Proof. Consider the function v £ (x) = u e e^ x '~ Pe ^ x ''' £ , where p £ = 2d(x) — (3d 2 {x)/e and /3 is a 
positive parameter to be chosen later. Since C £ u £ — \ £ u £ = we have 

= C E (e^~^ x ^ +Pe ^^ £ v £ ) — X e e^^ +Pe ^^ e v £ = e^^ +Pe ^^ e C e v e , 

the operator C £ being given by 

~ (fin f) f)u 

C £ v = ^{x,x/e-)^^ + e(b 3 (x,x/e a ) - 2a^ — (<j ) - Pe))^ + 

with 

c £ (x) = c(x,x/e a ) -X £ + H(V((f) - p £ ),x } x/e a ) - ea lJ (x,x/e a ) ^ — (0(x) - p 6 (z)). 

The coefficient c e (x) depends on the parameter /3 through the derivatives of p £ , and we can 
choose so large that c £ < in Q \ £l B /p for every e > 0, where f2 e /^ = {a; G f2; d(x) > e/fl}. 
Indeed, |Vp e | < 4 when d(x) < e//3 while A e > —A (cf. Lemma [7J and 



a iJ '(x, x /e ot )-^ r < --/3a ij (x, z/O a a + C ^ ~2m/3/e + C in Q \ n e/j3 , 



where C is independent of e and /3. Thus there is /3 > such that c £ (x) < in Q \ Q. e /p. With 
this choice of /3 we get by the maximum principle applied to the equation C £ v £ = 0, that v £ < 
max{v £ (x); x G dVL £ /p} in fl \ fl e//9 , or u £ e^ x)/e < e {p ^ x) ~ £/ ^ /£ max{u £ (x)e^ x)/e ; x G dQ e/p } 
when c?(x) < e/(3. Since p e (a;) < e/(3 when <i(x) < setting re = 1/(3 we get the desired 
bound. □ 

Corollary 10. There is a constant re > such that every maximum point x £ of u £ satisfies 
d(x e ) > KE. 

Proof. We simply apply Lemma [9] with = 0. □ 
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Lemma [S] and Corollary [TU] imply that we can extract a subsequence of functions W £ con- 
verging uniformly on every compact in Q to a limit W(x) G C(Q). Moreover, by Lemma [8] the 
function W is Lipschitz continuous and max{|W e (:r) — W(a;)|; i6 0, d(x) > ks} — > for every 

K > 0. 

Corollary 11. Let converge (along a subsequence) to a function W uniformly on every 
compact in Q. Then, for every G C 2 (Q), we have 

limmin{W / e (x) — <f>(x)\ x G Q} = min{H^(a;) — (f>(x); x G Q}. 

E— S-0 

Proof. We know that max{|W e (x) — iGll, <i(x) > ne} — > as e — >■ for every k > 0. 

On the other hand, by Lemma [U min{W / e (a;) — <f)(x);x G Q} is attained at a point located on 
the distance at least k, e from the boundary dQ for some k q > independent of e. This implies 
that liminf £ _ 5 .o mm {W / £(x) —<f>(x); x G fi} > min{W / (x) — a; G f2}. The opposite inequality 
limsup^o min{W / £ (x) — (p(x);x G Q} < min{W / (x) — (/>(x);x G fi} is an easy consequence of 
the uniform convergence of W £ to W on compacts in Q. □ 

5 Vanishing viscosity limit 

This section is devoted to the proof of Theorem [TJ According to the results of the previous 
section we can assume that, up to a subsequence, 

K -> A (5.1) 

and 

m&x{\W £ (x) - W(x)\; iffl, d(x) > Ke} -> > 0, (5.2) 

where W is some Lipschitz continuous function. We are going to show that the pair A and W 
is a solution of problem (jl.3p . (II. 4p . 

We follow the same scheme for a > 1, a = 1 and a < 1. 

For an arbitrary function <p G C 2 (f2), let W — (f) attain strict minimum at a point Xq G Q. 
We construct test functions of the form 

<f>e = 4>e-\x- X \ 2 , 

where e — > <p uniformly in Q. Moreover, we will chose functions <p £ with uniformly bounded 
first derivatives and such that 

-^(x £ , l£ / e tt )^(x e ) + i/(V^^),x E ,x E /£ a ) ^F(V0(x o ),x o ) (5.3) 

for every sequence of points x £ G Q such that x £ — > xq. The existence of such functions <p £ will 
be justified later on. 
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Let us now show that minima points of functions W £ — (f) £ converge to Xq. Indeed, 

| min(W £ - (f> £ ) - mm(W e — </>+\x — x \ 2 )\ -> 

since (p £ (x) — > (f)(x) — \x — Xq\ 2 uniformly in Q. If x £ is a minimum point of W £ — <p £ then in 
view of Corollary [11] we have 

min(W £ - (f) £ ) = W £ (x £ ) - <p(x £ ) + \x e - x \ 2 + 6(1) > min(W - 0) + \x e - x \ 2 + 6(1). 

On the other hand, choosing a sequence of points x £ G Vt such that d(x e ) = \x £ — xq\ = e we 
have 

min(W £ - </> e ) < W £ {x £ ) - e (x e ) -> W(zo) - <f>{x ) = min(I^ - 0), 
therefore x £ — > Xq. 

Now, if x £ is a minimum point of (W £ — <p £ ) then clearly x £ G Q, hence V<p £ (x £ ) = VW £ (x £ ) 
and a ^{x £ ,x £ /e a )^-{x £ ) > a^(x £ ,x £ /e a )^-(x £ ). We use (JO} to get 



2, 



ff(V^(i e ),a: e) j; e /e a ) - A £ > ea iJ '(x e> a; £ /0^— ^(x £ ), 



or 



- ea i i(x £ ,x £ /e a )-^^(x £ ) + H(V4> E (x £ ),x e ,x £ /e a ) - \ £ + 0{e + \x £ - x \) > 0. (5.4) 

Next we use (15. 3p to pass to the limit in (15. 4p as e — > that leads to the desired inequality 

#(V0(x o ),xo) > A. (5.5) 

If W — cj) attains strict maximum at a point x £ ^> we argue similarly. We construct test 
functions of the form £ = £ with £ satisfying (15. 3p and such that (j) £ — > <fi uniformly in Q. 
Note that this time Xo G" dQ and therefore we can always chose a sequence of local maxima 
points x £ of W £ — (f) £ converging to xq. Then using the same arguments as in the proof of ( 15. 5ft 
we derive 

H(V</>(x ),x ) > A. (5.6) 

Thus W(x) is a viscosity solution of ( II. 3p . ( II. 4p . 

It remains to construct functions <ft £ that have uniformly bounded first derivatives, satisfy 
( 15. 3 p and converge to uniformly in Q. 

Case a > 1. We set 

£ (x) = 0(x) +e 2Q - 1 fl(x/e a ), 
where 9{y) is a y-per iodic solution of 

- a**(x , y)-QjpQjp = H(V<f)(x ), x Q ) - H(V<f>(x ), x , y). (5.7) 
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Thanks to (12. 3p such a solution does exist, for ( 12. 3 p is nothing but the solvability condition for 
(15. 7p . Moreover, since the coefficients and the right hand side in (15. 7p are Lipschitz continuous, 
9 £ C 2,1 (see, e.g., |20]). Therefore if x £ — > Xq as e — > 0, then we have 

- ea«(x e , x £ /e a )^^- + # (V0 £ (x £ ), x £ , x e /e a ) = -a*(x e , Xe /e a )^^- + 0(e) 

+ H(V<P(x £ ) + 0(e a ~ 1 ),x £ ,x £ /e a ) = - a ^(x ,x £ /e a )^^- 

+ H(V(f)(x ) } x 0} x £ /e a ) +0(\x-x £ \+e + e a - 1 ) = H(V<f>(x ), x ) + o(l). 

Case a = 1. Set 4> £ (x) = <fi(x) + e9(x/e), where 9(y) = — log fl(y) and $(y) is the unique (up 
to multiplication by a positive constant) Y — periodic positive solution of 

ai ^ X(hV) Wdy- + V{v) W- + t{v))q& = (5 ' 8) 

where p = V(f>(x ), V(y) = b>(x , y) - 2a lJ (x , y)pu c(y) = a lJ (x ,y)piPj -b 3 (x ,y)pj + c(x ,y). 
By a standard elliptic regularity result we have 9 £ C 2,1 (see |20j). 
Simple calculations show that 9(y) satisfies 

- a^(x ,y)-^- + H(V(j)(x ) + V9 } x ,y) = H (V0(x o ), x ). (5.9) 

Then we easily conclude that 

- ea*(x g , x e /e) d 2f x e) + g(V0 £ (x £ ), x E , x e /e) = -a«(x e , x e /e)^^- + 0(e) 

+ H(V<f>(x E )+V9(x E /e),x E ,x E /e) = -a* fro, x e /e) 
+ H(V(f)(xo) + V9(x £ /e),x ,x £ /e) + 0(\x £ - x \ + e) = H(V<f>(x ),x ) + o(l), 
as soon as x £ — > Xo when e — )■ 0. 

Case a < 1. Set £ (x) = <p(x) + e a 9 £ (x/e a ) where 9 £ is a F-periodic solution of the equation 

- eW (x , y)^^ + # (/> + V0 e (y), z , y) = # e (p, ar„) with p = V0(x o ). (5.10) 

Such a solution exists if H e (p, Xq) coincides the first eigenvalue \x e (eigenvalue with the maximal 
real part) of the spectral problem 

e 2 ^a^(x , y)^: + e l - a V{y)-^- + S(l/)0 e = /^ £ , 
i? £ is ^-periodic, 
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where c are as in ( 15. 8p . According to the Krein-Rutman theorem fi £ is a real and simple 
eigenvalue and the corresponding eigenfunction d £ can be chosen positive. Then a solution of 
(I5.10P is given by 9 £ = —e x ~ a log-$ e . We invoke now the results obtained in Section [31 

H £ (p,x ) ^H(p,x ) = H(V(j)(x ),x ) (5.11) 

(where the limit H(p,xo) is described in (12. 5p ). 

||9 2 C/%%||l- <Cje l - a (5.12) 

This allows us to show ( 15.31) similarly to other cases considered above, 

- ea«(x e , x £ /e a ) ^* e) + #(V0 £ (a; £ ), x £ , x £ /e a ) = -e^V^, WO ^^ + £00 

+ H(V<f)(x £ ) + V6 £ (x £ /e a ), x £ , x £ /e a ) = -eW (ar , ^/e a )^^^ 
+ H(V(j)(x ) + Ve £ (x £ /e a ),x , x £ /e a ) + 0(|x - x e | + e) = #(V0(x o ), z ) + o(l). 
Theorem [1] is completely proved. □ 



6 Lower bound for eigenvalues via blow up analysis 

From now on we consider in details the special case when a = 1 and c(x, y) = 0. Under the 
assumption (12. 17ft the eigenvalues A e of ( 11. ip converge to zero as e — > 0. We are interested 
in the more precise (up to the order e) asymptotics for X £ . We resolve this question by local 
analysis near points £ of the Aubry set Ajj of the effective Hamiltonian. 
Fix a point £ G Ajf. 

Applying the maximum principle we see that X £ < 0. On the other hand, it is well known 
that the eigenvalue X £ is given by 



X £ = inf (sup -4~; G C 2 iSl) n C(n), > in fi, = on dtt\. (6.1) 



Therefore, for every given 5 > 0, we have 



A e > ecr e , 



where a £ < is the first eigenvalue of the equation C £ v £ — 5\x — £| 2 f e = £<5" £ t> £ or 

ea lJ (x,x/e) ^ e + V{x,x/e)^- - — - - t> £ = a £ t; e in ft (6.2) 

with the Dirichlet condition w e = on <9f2. We assume hereafter that the first eigenfunction v £ 
is normalized by f e (£) = 1. 
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Let us transform (16.21) to a form more convenient for the analysis. First, after changing 
variables z = (x — £)/y/£ an d setting w £ (z) = v £ (^ + \fsz) equation (16. 2p becomes 

air i r\ d2w z , b u/ei^ z ^ z /\^) dw £ , . ^ 

a L/e\V £Z i Z /v £ )-Q^-. + 7y= ~ 5 N W £ = a eW e IB (ft - (6-3) 



Here (and below) the subscript "£,£/£" denotes the shift (translation) by £ in x and by £/e in 
y, i.e., for instance, a l ^, £ (x,y) = a l ^, £ (£+x,£,/e+y). Next multiply (|Q]| by 9*^ e (y/ez, z/y/e), 
9*(x,y) being given by ( 12. 16ft . to find after simple rearrengements 

d fa* { r i r\ ij < r i r\ dw A , s u/e(,V^ z i z /V^} dw e 



+ ( &J(v ^ + ° + ^W)^ = (*« + *|*| 2 )*k, 6 (V«, VVSR (6-4) 
where S 3 ^, £ (x,y) is obtained by shifts (as described above) from 

S\x, y) = V{x, y)9*(x, y) - A (a^'(x, (x, y)) - F(x), 

and are uniformly bounded functions. Since 9* solves ( 12 . 1 6j) . the F-periodic vector field 
S(x, y) = {S 1 (x, y), . . . , S (x, yj) is divergence free, for every fixed x, and (due to the definition 
of b) this field has zero mean over the period. Therefore we can find the representation 

8 

S J (x,y) = -g—T lJ (x,y) with F-periodic in y skew-symmetric T lJ (x,y) (T 13 = —T ]l ). 

Moreover, functions T %3 are (can be chosen) continuous with bounded derivatives dT 13 /dxu- 
We can thus rewrite (16.41) as 

A (g* /e (Viz, z/yfi ^) + ( *' ( ^ + S + V?h>(zj) |g = (a £ + 5|,| 2 )^ /£ (^,^)- £ , 

_(6.5) 

where q^,(x,y) = + x,£/e + y), q tJ (x,y) = 9*{x 1 y)a l] (x,y) + T tJ (x,y), and h{ are 

uniformly bounded functions. Note that on every fixed compact we have 

uniformly in z as e — > 0. 

Lemma 12. //&(£) = /or some £ E ft £/ien £/ie /irst eigenvalue X £ of the operator {2. 14 ) 
satisfies the bound —Ae < X e < with some A > independent of e. 

Proof. We know that A e < and in the proof of the lower bound we assume first £ G ft. Then 
(16. 5 P holds in B 2 = {z; \z\ < 2} for sufficiently small e. Let us write (16. 5p in the operator form 
Cs We — ^ £ 9^/ £ (y/ez, z/^)w £ and consider the parabolic equation for the operator £, £ aux ^ 

dt 
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subject to the initial condition w £ (0,z) = w e (z) and the boundary condition w E (t,z) = 
on (0, +00) x dQ. The solution w £ of this problem has the pointwise bound w £ (t,z) < 
exp(a £ (min#*)t) w E (z). This follows by the maximum principle applied to 



wJz) > 0. 



On the other hand, since the coefficients of the operator C £ are uniformly bounded and the 
uniform ellipticity bound 6^i E [y/ez, z/ y/e)a % ^,{y/ez, z/\/i)C*C" 7 ' > (min 6>*)m |£| 2 holds, by 
the Aronson estimate (see |3j) we have 

mm{w e (l, z); z G B{\ > M min{w e (0, z); z G B\) 

with M > independent of e, where B\ is the unit ball B\ — \z\ \z\ < 1}. This yields 

e ( mine *)^ minw £ > mm{w £ (l,z); z G B x } > Mmin{w e (0, z); z G B x ) = Mining, 

i.e. a £ > logM/min6>* =: -A. Thus a £ > -A and A e > -Ae. 

Finally, in the case £ G dfl we can repeat the above argument taking £ e G f2 in place of £, 
with |£ 6 - £| = dist(f e , dVl) = 2y/e. □ 

In the proof of Lemma fT2l we have got a uniform lower bound for a £ which (in conjunction 
with the obvious inequality o £ < 0) allows one to obtain uniform bounds for the norm of w £ 
in C°'"(K) (with (3 > depending only on bounds for coefficients in (16.5)) ) and if 1 (i^), for 
every compact K (see, e.g., |20|). Thus, up to extracting a subsequence, w £ — > w in Ci oc (IR w ) 
and a £ — > a. Moreover, using well established homogenization technique based on the div-curl 
Lemma we get that w solves 

Q ij ^- + B i ^z i ^-5\z\ 2 w = aw in R N , (6.6) 

OZiOZj OZj 



where = Qi l are some constant coefficients satisfying the ellipticity condition (actually, one 

1 d 2 H 

2 dpidpj 



can check that Q ij = |/#-(0,£)) & nd 



BH= B^(0) = ^(0- 

Since we assumed the normalization w £ (0) = 1, we see that w(z) is a nontrivial solution of 
(16. 6p . Moreover, if z £ is a maximum point of w £ (z) we get from (16. 3p \z £ \ 2 < —o £ jb therefore, 
thanks to Lemma [T2| \z £ \ < C. It follows that w(z) is a bounded positive solution of ( 16. 6p . 
A solution of (16. 6p can be constructed in the form w(z) = e~ v ^ ZiZ i with a symmetric positive 
definite matrix (^s)ij=TN' Indeed, consider the following matrix Riccati equation 

ATsQTs - BTs - T 5 B* -51 = 0, 
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where / denotes the unit matrix. It is well-known that there exists a unique positive definite 
solution T$ (since 5 > and Q is positive definite, see [28]). Then w(z) = e~ v ^ ZiZ i is a positive 
bounded solution of (16. 6p corresponding to the eigenvalue a = — 2tr(Qr,5). Next observe that 
by means of the gauge transformation w(z) = e~ r ^ z ^w(z) (r > 0) equation (16. 6ft is reduced to 

g^T— — + (B 13 + 4rQ lj ) Zi — + (Ar 2 Q i3 z iZl + 2rtrQ + 2rB ij z i z j - 5\z\ 2 )w = aw in R N . 

OZiOZj OZj 

For sufficiently small r > we have ((4r 2 Q tJ Z{Zj + 2rtrQ + IrB^ZiZj — 5|-z| 2 ) — > — oo and 
w(z) — > as \z\ — > oo. Then, according to [37], the eigenvalue a corresponding to such 
a positive eigenfunction w (vanishing as \z\ — > oo) is unique. Thus o = —2tT(QTs), and 
summarizing the above analysis we have liminf^o A e /e — —2tr(QTs). Finally note that T$ 
converges to the maximal positive semi-definite solution of the Bernoulli equation (see, e.g., 

ATQT -BY - TB* = 0, (6.7) 

as S — > +0. Calculations presented in Appendix ICl show that — 2tr(Qr) = cr(£) with <j(£) being 
the sum of negative real parts of the eigenvalues of —!?(£). Thus, after maximizing in £ G Ajj 

we get 

liminf \ £ /e >W= max{(x(£); £ G Aw}- (6.8) 

7 Upper bound for eigenvalues and selection of the additive 
eigenfunction 

In this section we derive an upper bound for eigenvalues which completes the proof of the 
asymptotic expansion formula ( 12. 18[) . Similarly to the previous section we make use of the blow 
up analysis near points of the Aubry set. We consider here only special (so-called significant) 
points of the Aubry set, where we can control the asymptotic behavior of rescaled eigenf unctions 
at infinity. We will show that only these special points matter to the leading term of eigenvalues 
and eigenfunctions. 

Due to Theorem (TJ up to extracting a subsequence, the functions W £ = — e\ogu £ converge 
uniformly on compacts to a viscosity solution W of H(VW(x), x) = in Q, H(VW(x), x) > 
on dfl. It follows that W has the representation W(x) = minjfi^x, £) + W(£,); £ G Ajj}. 

We will say that point £ G Ajj is significant if 

W(x) = dj[(x,l;) + W(£) in a neighborhood of £. 

Otherwise we call £ negligible. For every negligible point £ G A-jj there are sequences x k — > £ and 
£ fc G Ajj\ {£} such that c%(x fc , f ) + W{£) > d w {x k , £ fc ) + W(£ k ). Passing to the limit (possibly 
along a subsequence) and using the continuity of the distance function we get d-jj (£, £') = 
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- W(?) for some f G Ag, ? ^ f (we always have > ^(0 - Now let 

us introduce the (partial) order relation ^ on .4.^ by setting 

r d e = ^(0 - w{£). (7.i) 

This relation is clearly reflexive, its transitivity is a consequence of the triangle inequality 
dir(£, 6") < ^h(C) O + £") while the antisymmetry follows from the inequality Sjj(£, £') > 

held for all £, £' G Ajj with £ 7^ Then we see that every minimal element £ G .4# is a 
significant point. Since Ajj is finite there exists a minimal element, i.e. there is at least one 
significant point £ G 

Let us fix a significant point £ G ./4#. From now on we assume that u £ is normalized by 
w e (£) = 1, unless otherwise is specified; the W will also refer to the limit of scaled logarithmic 
transformations of u £ normalized in this way. Thanks to the upper and lower bounds for the 
eigenvalue \ e we have X £ /e — > a as e — > along a subsequence. Then we argue exactly as in 
the proof of the lower bound for A e . We consider rescaled eigenfunctions w £ (z) = w e (£ + y/sz) 
that are solutions of 

Up to extracting a further subsequence, they converge in C(K) and weakly in H l (K) (for every 
compact K) to a positive solution of 

Q ij + B ij z~ = aw mR N . (7.2) 

OZiOZj OZj 

Eigenvalue problem (17. 2p possesses, in general, many solutions even in the class of positive 
eigenfunctions w(z). We will uniquely identify a studying the asymptotic behavior of w(z) as 
\z\ — > 00. More precisely, we will show that 



w(z 



i e M|n*z| 2 v\n* u z\ 2 ig b ounc i ec i on £ Qr gome fx> an( j ever y u > 0, (7.3) 



where U s and U u denote spectral projectors on the invariant subspaces of the matrix B corre- 
sponding to the eigenvalues with positive and negative real parts (stable and unstable subspaces 
of the system Zi = —B^z/). This allows to use the following uniqueness result. 

Lemma 13. Let w(z) be a positive solution of 117. S\) satisfying A7.3\) . Then w(z) = Ce -r * ZiZj 
with C > 0, where F is the maximal positive semi-definite solution of W. 7\) . Moreover, we have 

a = -2tr(rg). 

Proof. First observe that w(z) = Ce~ T ZiZi satisfies ( 17. 3p . This follows from the relation T = 
n s rn* > 7n s Il* where the inequality holds for some 7 > in the sense of quadratic forms, see 
Proposition [T7] in Appendix ICl It is also clear that w{z) does solve (17. 2p with a = — 2tr(rQ). 
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To justify the uniqueness of a and w(z) we make use of a gauge transformation w(z) = 
e^ z >w(z), with a quadratic function <p(z) to be constructed later on, which leads to the equation 
of the form 

cP 1 in ~ c)i 11 ~ 
Q ij -^^ + B»z i — + C(z)w = aw in R N . (7.4) 

We will choose <p(z) so that C(z) — > — oo, w(z) — > as \z\ —> oo. Then, by [37J, there is a 
unique a such that (I7.4p has a positive solution w(^) vanishing as \z\ — >■ oo (w(z) is also unique 
up to multiplication by a positive constant). 

Construct <j)(z) by setting <ft = rA l J ZiZj — rA^ZiZj, with symmetric matrices A s and A u , to 
get B ij = B*j + ArQi\A l l - A h s ) and 

C{z) = Ar\A l u - A u s )Q lm {A^ - A^)z iZj 

'(B u {A l Z - A l i) + (A^ - A»)B l i)z iZj + 2tr(Q(A u - A s ))). 

Define A s and A u as solutions of the Lyapunov matrix equations 

BA 8 + A S B* = n s n:, BA U + A U B* = -U U U* U , (7.5) 

given by 

e Bt U s U;e BH dt, A u = e Bt Ii u W u e BH dt, (7.6) 

oo JO 

and choose sufficiently small tq > in such a way that the matrix 

Ar(A u - A S )Q(A U - A s ) + B(A U - A.) + (A u - A S )B* = Ar(A u - A S )Q(A U - A 8 ) - U S U* S - U U U* U 

is negative definite for < r < r . Then C(z) — > — oo as \z\ — > oo. It remains to see that if 
w(z) satisfies ( I7.3P then choosing small enough r > we have w(z) — > as \z\ — > oo. Here we 
have used the inequalities A s < 7iII s n* and A u > 72il u n* for some 71,72 > 0. □ 

So far we know that X £ /e — > o and w £ (z) = w e (£ + \fez) converge in uniformly on compacts 
to a positive solution of ( 17. 2p . In order to apply Lemma [TBI we need only to show (17 .3j) . To this 
end we first construct a quadratic function ^(x) satisfying 

H(V^(x),x) < -8\x - £| 2 in a neighborhood U(£) of £ (7.7) 

for some 5 > 0. 

Lemma 14. Let us set <p s (x) := A l Jx{Xj and <p u (x) := A^XiXj, where A s and A u are solutions 
of the Lyapunov matrix equation ( 17. 5\) given by Ii7.6]) . Then the function 

:= ii<t> s (x - - v<t>u{x - 
satisfies ( 17. 7\) for some 5 > 0, provided that < fi, v < r and r > is sufficiently small. 
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Proof. We have, as x — > £, 

_ _ BH d§ u 

H(V^(x),x) < H(0,x) + —(0,x)-^(x) + C\V^(x)\ 2 

= -fa - 6)^(0 ^(*) + C|V$;(x)| 2 + o(\x- £| 2 ) 
< -2( Xi - - vA»){ Xl - 6) 

+ d^in:^ - or + ^ra* - on + * - o 2 )- (7.8) 

Note that B{^A S — vA u ) + (fj,A s — vA u )B* = fiU s Il* + z/ILjIT*, therefore the first term in the 
right hand side of f!7.8p can be written as /i|n*(x — 0| 2 + ^|n*(a; — OP- Thus f!7.7p does hold 
if < \i < l/Cx and < v < 1/d. □ 

Next we prove 



Lemma 15. If <& v {x) and //, v are as in Lemma\n\ then W(x) > 3>„(x) in U'(£) \ where 
U'(£) C £/(0 i s a neighborhood of £. 

Proof. Since £ is a significant point, we have W(x) = djj(x,£) in U(£). Due to the represen- 
tation formula ( 12. lip <%(x,0 = lim^oo f Q L(f] k ,r] k ) dr for a sequence t k > and absolutely- 
continuous curves i] k : [0,t k ] — > Q satisfying the initial and terminal conditions rj k (0) = £, 
r] k (t k ) = x. We claim that there is a neighborhood U'(£) C £7(0 such that Vx G E/'(£) we have 
r} k {r) G £7(0 £ [0,^ fc ] when /c is sufficiently large. Indeed, otherwise there are sequences of 
points x k — > £ and curves f) k (t), rj k (0) = £, r] k (t k ) = x k that exit t/(0 a ^ a time t = r k and 
lim^oo J L(r] k , rj k ) dr = 0. Let us set y fe := 7? fc (r fc ) G dU'(£), then hm^oo Sjj(y k , = and 
after extracting a subsequence y k — > y G dU'(£) we obtain S-jj(y,£) = 0. Therefore y G A^-. 
We can repeat this reasoning to find y G *4.#n dU"(£) for every open subset £/"(£) of £/(0 
containing point £. Thus £ cannot be isolated point of ^4-^-, contradicting ( 12. 17ft . 
Now using (17.7)) we get, for every a; G £/'(0 

*£(x) = f V^(i7*)-i)*dT= f (V^( V h )-rj k -H(V^(v k ),V k ))dr 
Jo Jo 

+ f H(V^( V k ), V k )dr< f Ltf,rf)dT, 
Jo Jo 

when k is sufficiently large. It follows that $^ < W in U'(£). On the other hand if $^ = W 

at a point Xo G t/'(0 then Xo is a local minimum of W — $> v and if (V$^(xo), Xo) > since W 

is a viscosity solution of /f(VW(x), x) = in f2. Therefore xq = £ by (17. 7ft . i.e. <£>^ < in 

£7'(0 \ {£} an d by choosing, if necessary, a smaller U'(£) we are done. □ 

The following is the crucial step in establishing (17. 3p . We construct a test function ^ £ (x) 
of the form \l/ e (x) = $^(x) — e6 £ (x,x/e) which satisfies 

- ea ij (x, x/e) + (W £ (x), x, x/e) < H(V$"(x),x) + Ce mU'(£). (7.9) 

OXiOXj 
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Assuming first that the solution i9(p,x,y) of ( 12. 4p . normalized by J Y , d(p,x,y)dy = 1, is suf- 
ficiently smooth, we set 9 £ (x,y) = #(V$^(x), x, ?/), where 9(p,x,y) = \og{}(p,x,y). Then, 
since 

-a %3 {x,y) Qy' d y. +H[p + V y 6(p,x,y),x,y) = H(p,x), 

one easily shows (17. 9p . Note that in this case 8 £ (x,y) is independent of e. In the general 
case, thanks to (^-regularity of the coefficients a lJ (x,y) and V(x, y), all the first and second 
order partial derivatives of i9(p, x, y) exist and continuous on R N x Q x ~K N , except (possibly) 
d 2 $(p,x,y)/dxidxj. To obtain sufficient regularity of 8 e (x,y) we set 



6 e {x,y) = f <f E (x-x')9(V$»(x),x',y)dx', 



where (p £ (x) = e <p(x/e), with <f(x) being a smooth compactly supported nonnegative function 
and J <f(x) dx — 1. Then we have 



and \V y e(V^{x),x,x/£) - V (6 £ (x,x/e)) \ < C. This eventually leads to flZjfl . 
It follows from (EHJ) and ([72]) that 



d 2 \fj 

-ea l \x,x/e) gx e + H(VV e (x), x, x/e) < -5\x - £| 2 + Ce in [/'(£)■ (7.10) 

Consider now the function W £ — ^> £ . By Lemma [151 we have W £ > \l/ £ on dU'(£) for sufficiently 
small e, therefore either W £ > ty E in £/'(£) or W £ — \fr e attains its negative minimum in £/'(£) 
at a point x £ . In the latter case we have WW £ (x £ ) = W^ £ (x £ ) and a 13 (x £ , x £ /e)^^-(x £ ) > 
a t3 (x £ ,x £ /e)^^-(x £ ) 1 this yields 

d 2 W £ 

X £ = -ea lJ (x £ ,x £ /e) — ^ (x e ) + # (VW e (a; e ),ar e ,ar e /e) 

d 2 ^ 

< -ea l3 (x £ ,x £ /e)——^(x £ ) + H(W £ (x £ ),x £ ,x £ /e) < -5\x £ - £| 2 + Ce. 

OXiOXj 

Thus either W £ > ^ £ in [/'(£) or W £ > ^ £ + W £ (a; e ) - ^ £ (a; £ ) in [/'(£) and x £ satisfies |a; e -£| < 
C-^/i. Both cases lead to the bound W £ (x) > 3>„(x) + W £ (x £ ) — (3e, where x £ is either £ or x £ 
(recall that it e is normalized by w £ (£) = 1, i.e. W e (£) = 0). Then setting z = (x — £)/y/e we get 

w e (z) < Cw £ (z £ )e-^ + ^ in ([/(£) - 

where z £ = (x £ — and hence |z e | < C. Observe that since z £ stay in a fixed compact as 

e — > 0, then u> £ (.2 £ ) < C and in the limit we therefore obtain 



w(z) < Ce-^ s{z)+u<t>u{z) in R N . 
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It remains to note that <f) s (z) > 7 3 |I1*2;| 2 and <p u (z) < 7 4 |n*z| 2 for some 73,74 > 0. Hence w(z) 
does satisfy (17. 3p . and by Lemma [T3l we have a = —2tr(TQ) = cr(£). Thus 

limsup X £ /e = <r(£) for every significant point £. (7-11) 

£-S>0 

Inequalities (16. 8p and (17. lip prove formula (I2.18p . Moreover they imply the uniqueness of 
the limiting additive eigenfunction W(x), provided that the maximum in (I2.18P is attained 
at exactly one £ = £. Indeed, we know that, up to extracting a subsequence, functions W £ 
converge uniformly (on compacts in Q) to an additive eigenfunction W(x); here W £ = —elogu £ 
and u £ are referred to the eigenf unctions normalized by (12. 2p . By (17. lip the unique significant 
point (associated to the chosen subsequence) is £. Therefore £ is the only minimal element in 
A-jj with respect to the order relation -< defined in (17. ip ; hence it is the least element of Ajj-, i.e. 
£ ^ £ for every £ G Ajj. This means that W(£) = W(£,) + (%(£, £) V£ G Ajj, and consequently 
W(x) = djj(x,£) + W / (£). Finally, by Corollary (TT] we have W / (£) = 0, and Theorem [2] is now 
completely proved. □ 



8 Other scalings 

Theorem [2] can be generalized to the case of 5 a -scaling, a > 1, of the fast variable in (II. ip 
(with c(x,y) = 0). In this case the effective drift is still defined by formula (I2.15P where 6* is 
now the F-periodic solution of the equation d ^ dy (a y (x, y)9*(x, y)) = normalized by J y 8* dy. 
However, the discontinuous dependence of the effective drift on the parameter a > 1 (at a — 1) 
might lead to a significant shift of the concentration set of the eigenfunction u e from points £ 
where 6(£) = 0, if a > 1 sufficiently close to 1. 

We outline main changes to be made in order to adapt the arguments of Sections [6] and [7] 
to the case a > 1. First of all let us introduce the approximate Hamiltonian H e (p,x) as the 
(additive) eigenvalue corresponding to a ^-periodic eigenfunction of 

- a ij (x, y) d26 ^ v) + H(p + e^Vyd^p, x, y),x, y) = H £ (p, x), (8.1) 

and the approximate drift b £ (x) by 

bi(x) = -^(0,x). 

The eigenvalue H £ is unique and 6 £ is unique up to an additive constant, moreover 9 £ can be 
found as the scaled logarithmic transformation 6 £ = — 2 (i-i) log^e of a positive F-periodic 
eigenfunction of the problem 

,92 o a a 

s 2 ^a ij (x,y)-^- + e 1 - a (V(x,y) - 2a* (x, y)Pi)^ + H(p, x, y)$ £ = H £ (p,x)# £ . 



24 



Similarly to the case a = 1 the drift b e (x) can be equivalently defined by 

b e (x) = J b(x,y)9*(x,y) dy, 
via the F-periodic solution 6* of the equation 

It is clear that b £ — > b in C 1 topology, therefore if b has, say n, zeros in Q, where the 
hyperbolicity condition (for the ODE x = —b(x)) is satisfied, then b e has exactly n zeros at the 
distance at most 0{e OL ^ 1 ) from the corresponding zeros of b. Then, to show the lower bound for 
eigenvalues one follows the lines of Section |6] with a zero £ £ of b e in place of the corresponding 
zero £ of b, and 8* in place of 9*. Note that although £ e — > £ as e — > 0, the distance between this 
two points might be of order e"" 1 , so that in the local scale y/e this distance tends to infinity. 
Nevertheless, up to the shift from £ to £ £ the local analysis is exactly the same as in Section |6j 
Let us emphasize that for a 6 (1,3/2) the statement of of Lemma fT2l remains valid only if at 
least one of zeros of b is an interior point of Q. 

The argument of Section [7] can also be adapted to the case a > 1. As in the proof of the 
lower bound one finds equation ( 17.2ft for the limit of rescaled functions w £ (z) = w e (£ e + \/sz), 
while the construction of the functions $^ and \l/ e is to be modified. One can linearize the drift 
b £ at £ e and construct the quadratic function (which now depends on e) following Section [7] 
with B l J = §|:(£e) in place of B^; also, in the construction of the function \I/ e one makes use 
of the eigenfunction 6 e (cf. (jSTTj) ) and sets * £ (x) = $^(x) + £ 2a_1 £ (V$^(a;), x, x/e a ). Details 
are left to the reader. 

Finally note that the case a < 1 remains completely open. 

9 Example 

Here we consider an example of an operator of the form (12. 14ft for which conditions ( 12.17P are 
fulfilled. Let x y (t) be a solution of the ODE x y = —b(x v ), x y (0) = y. We assume that 

• The vector field b(x) has exactly three zeros £ , £ 2 , £ 3 in Q. All of them are interior points 

of a 

• £ x and £ 3 are stable hyperbolic points, that is eigenvalues of ( — ||-(£ 1 )) i j-yn an( ^ 
( — lx"^ 3 ))i -Tlv ^ ave negative real parts; £ 2 is a hyperbolic point and cr(£ 2 ) > 
max^Vte 3 )}. 

• The ODE x = —b(x) does not have a solution with \im t ^ +QO x(t) = lim<^_ 00 x(t) = £ 2 . 

• For every y e U\ ULi{^}, either lim x y (t) = £ 2 , or inf H < ; x y {t) G H} > -oo. 

•> t— >— oo 
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One can easily check that under these assumptions the Aubry set Ajj coincides with lJj=i{£ J }- 
Hence, by Theorem [2j W(x) = c%(x,£ 2 ) and A e = ecr(£ 2 ) + o{e). 

It is interesting to trace in this example the possible structure of the set Z = {x G 
Q ; W(x) = 0}. It depends on whether there are trajectories of the equation x = —b(x) 
going from £ x or £ 3 to £ 2 , or not. 

Let Z 1 be the set of all points y G fl such that lim x y (t) = and lim x y (t) = £ 2 , and 

t— 5>+oo t— >— oo 

let Z 3 the set of all points y G Q such that lim x y (t) = £ 3 and lim x v (t) = £ 2 . It follows 

t— >+oo t— >— oo 

from (l2~XTj) that Z = {£ 2 } U^uF. 



10 Appendices 

A Uniqueness of additive eigenfunction 

The following simple result is a uniqueness criterion for problem (11.3j) - (ll.4j) . 

Proposition 16. Let A = so that $1.3\) - (L4\ ) has a solution W . Then W is unique (up to 
an additive constant) if and only if S H -_ x (x,y) = for all x,y G Ajj_ x , where Sjj_ x (x,y) = 
d H-\(x,y) + d-H_ x (y,x). 

Proof. If Sjj_ x (x, y) = then W(x) — W(y) = d H -_ x (x, y), since Wx, y we have W(x) —W(y) < 
djj_ x (x ) y). In particular, if Sjj_ x (x,y) = Wx,y G Ajj_ x then pick £ G Ajj_ x to get 
W(x) = d-jj_ x (x,() + W{£) on Ajj_ x . Thus W(x) = %_ A (x,£) + in (according to 

the representation formula ( 12.131) ). i.e. W is unique up to an additive constant. 

If there are two points £, £' G Ajj__ x such that Sij\_ A (£, £') > 0, then Wo(^) = djj_ x (x,!;) 
and W^(x) = %_ A (z,£') - c%_ A (£,£') are two solutions of (Ojl -lO )) and = W„(f) = W^g), 
while W (O - Wx(?) = %_ A (£, £') > 0. □ 



B Aubry set for small perturbations of a gradient field 

We outline here the proof of the claim stated in Remark[3j Consider a vector field b(x, y) which 
is a C 1 -small perturbation of VP (a;), i.e. \\b(x, y) — VP(x)||ci = 5 and 5 is sufficiently small. 
Let us show that the Aubry set Ajj of the Hamiltonian H(p, x) given by (12.41) (with c(x, y) = 0) 
is exactly the set of zeros of b(x) in Q, provided that 5 is sufficiently small and P G C 2 (Q) is 
as in Remark [3j 

Without loss of generality we can assume that H(p,x) = ^2,p 2 —V(x)pi, since the Aubry 
set of this Hamiltonian coincides with that of the effective Hamiltonian given by ( 12. 4ft . Let us 
first find the Aubry set Aho of the Hamiltonian H°(p,x) = J^pf — p^ ■ We calculate the 
corresponding Lagrangian L°(v,x) = \\v + VP(x)| 2 and use criterion (I2.12p . Let £ G Aho, 
then there exist a sequence of absolutely continuous curves rj n : [0, t n ] — > f2, r) n (0) = r] n (t n ) = £, 
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such that t n — > oo and lim^oo J Q * \fj n + VP(^ n )| 2 dr = 0. This yields 

0=lim / (\r] n \ 2 + 2\/P s (r] n )-r] n + \VP(v n )\ 2 )dT= Mm / (|?f| 2 + |VP(r/ n )| 2 )dr. 

Therefore, rj n (t) — > £ uniformly on every fixed interval [0, T\. It follows that £ belongs to the set 
K = {x G f2; VP(x) = 0}. Clearly, we also have K C Ajjo- Now note that the effective drift 
b(x) given by (I2.15p . can be written as b(x) = VP(x) + bs(x) with C 1 -small b$(x), \\bs\lc 1 — 0(5) 
as 5 — > 0. Thanks to the assumption on critical points of P(x), zeros of b(x) are isolated and 
are close to K when 5 is sufficiently small. Moreover, if u is a small neighborhood of £ G if 
then b(x) vanishes at exactly one point and |£ — = 0(5). Therefore, we can define 

a C 2 function P$ such that |VP,5(x)| > in Q \ K$, where K$ is the set of zeros of b(x), and 
\b(x) — VPg(x)\ = gs(x)\ VPj(x)| with max ;EG Q5f ( 5(x) = 0(5) as 5 — >■ 0. This yields the following 
bound (for small 5) 

\v + b{x)\ 2 >\\v\ 2 + 2WP&(x)-v + V s (x) } \/veR N , x G tt, 

where V§ > in Q \ K$. Then, arguing as above we see that Ajj = K$. Moreover, every £ G K$ 
is a hyperbolic fixed point of the ODE x = —b(x), as 5 is sufficiently small. 

C Properties of solutions of Bernoulli matrix equation 

We provide here some results on Bernoulli equation (16. 7p . used in Sections [6] and [71 Recall that 
the matrix Q in (16. 7\i is positive definite, U s and U u denote spectral projectors on the invariant 
subspaces of the matrix B corresponding to eigenvalues with positive and negative real parts. 

Proposition 17. The maximal positive semi- definite solution V of W. 7)) possesses the following 
properties: (i)T = n s ITl* ; (ii) V > 7ll s n* (in the sense of quadratic forms) for some 7 > ; 
(Hi) 2tr(Qr) = tr(PII s ), i.e. 2tr(Qr) is the sum of positive real parts of eigenvalues of B. 

Proof. It follows from (E2D that X = H U TU* U satisfies 

4n u rgrn; - (bu u )x - x(Bii u y = o. (c.i) 

Consider the symmetric solution of ( IC.lj) given by 

POO 

X= Y(t)dt, (C.2) 
Jo 

where Y(t) = -4e Bn ^n u rQrn;e( Bn ")*' (note that Y(t) = (BU u )Y(t) + Y(t)(BU u )* and 
Y(t) as t ->• +00, therefore integrating we get [BU U )X+X{BU U )* = -Y(0) = 4n„rQrn;, 
i.e. X does solve (IC.ip ). We claim that X = X. Otherwise Z := X — X is a nonzero solution 
of equation (BU U )Z + Z(BU U )* = and Z = Yl u ZYl^. Then Z(t) = Z is a stationary solution 
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of the differential equation Zit) = (BH u )Z(t) + Z(t)(BH u )* . The latter equation has the solu- 
tion Z{t) = e BIIll *n u Zn*e^' Bnt '' ) ** which vanishes as t — > +00 and satisfies the initial condition 
Z(0) = Z. Thus Z = 0, i.e. X = X. On the other hand it follows from (l02|) that X < while 



X > 0, this yields X = X = 0. Since F is positive semi-definite we also have H U T = ITI* = 
and the calculation T = (H u + n s )r(II u + II S )* = n s ITI* shows (i). As a bi-product we also 
have established that T is the maximal positive semi-definite solution of 

ATQY - (BTv s )T - Y{BU 8 )* = 0. (C.3) 



Indeed, assuming that T is another positive semi-definite solution of (I C 3 [) we get n u rQrn* = 0. 
This yields fn* = so that f = n s fn*, therefore BY = 5(II s ) 2 fn* = (BU s )f and f thus 
solves ( 16. 7ft . 

To show (ii) and (iii) consider the maximal positive definite solution T$ of 

4f sQT s - (BU S + 6I)f s - f 5 (BU s + 51)* = (C.4) 

for S > 0. The existence of the unique positive definite solution follows from the fact that Tj 1 
is the unique solution of the Lyapunov matrix equation 

4Q - fj\BU s + SI) - (BU S + Siyfj 1 = (C.5) 

given by 



f -1 = 4 I e (BU a+ Sirt Qe (BU 3+ 6I)t dt 
J —00 



It is known (see [28]) that converges to the (maximal positive semi-definite) solution V of 
( 1C.3[) as 5 — > +0. This allows to establish (iii) easily, 



2tr(Qr) = 2 lim 2tr(Qf s ) = - lim tr (ttHBUs + 6I)T S + (BU S + 51)*) = tr( BU a ). 

Finally, if we assume that (ii) is false, then there is 77 G W N such that T?7 = while 
7^ 0. Thanks to (i) the equality Trj = implies that m*^ = 0. On the 
other hand, r ^lim^ + o f J 1 H s H*r] S j = n s II*?7, where the limit lim^^+o f J 1 H s H*r] exists, for 
^BJis+siytQ^BYis+si^Yiji*^ decays exponentially fast as t — > —00, uniformly in 5 > 0. Ac- 
cording the Fredholm alternative n s II*?7 and 11*77 m ust be orthogonal, yielding | IT*77 1 = 0. We 
obtained a contradiction showing that (ii) does hold. □ 
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